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The theoretical Green-Kubo relation for viscosity is tested using experimentally obtained data. 
In a dusty plasma experiment, micron-size dust particles are introduced into a partially-ionized 
argon plasma, where they become negatively charged. They are electrically levitated to form a 
single-layer Wigner crystal, which is subsequently melted using laser heating. In the liquid phase, 
these dust particles experience interparticle electric repulsion, laser heating, and friction from the 
ambient neutral argon gas, and they can be considered to be in a nonequilibrium steady state. Direct 
measurements of the positions and velocities of individual dust particles are then used to obtain 
a time series for an off-diagonal element of the stress tensor and its time autocorrelation function. 
This calculation also requires the interparticle potential, which was not measured experimentally, 
but was obtained using a Debye-Hiickel-type model with experimentally determined parameters. 
Integrating the autocorrelation function over time yields the viscosity for shearing motion amongst 
dust particles. The viscosity so obtained is found to agree with results from a previous experiment 
using a hydrodynamical Navier-Stokes equation. This comparison serves as a test of the Green-Kubo 
relation for viscosity. Our result is also compared to the predictions of several simulations. 

PACS numbers: 52.27.Lw, 52.27.Gr, 66.20.-d, 83.85.Jn, 05.60.Cd 



I. INTRODUCTION 

There are various two-dimensional (2D) physical sys- 
tems that allow direct observation of individual parti- 
cle dynamics. These systems include colloidal suspen- 
sions [l[, granular materials 0, electrons on a liquid he- 
Hum surface in a Wigner lattice [1] , ions confined magnet- 
ically in a Penning trap 0], and single-layer dusty plas- 
mas In these systems, the relevant particles collide 
with their neighbors frequently, so that momentum and 
energy are transported from one place to another. (In 
all these 2D physical systems, motion is not purely 2D, 
but usually includes some limited out-of-plane motion, 
so that the systems are often described as quasi-2D.) 

Shear viscosity, 77, is a transport coefficient that char- 
acterizes the momentum flux perpendicular to a velocity 
gradient. Sustaining the velocity gradient requires the 
application of a shear stress, which corresponds to an 
ofF-diagonal element of a stress tensor. The hydrody- 
namical definition of shear viscosity is the ratio of this 
ofF-diagonal element and the velocity gradient As a 
measure of dissipation, viscosity is useful, for example, 
in describing the damping of shear waves 

The Green-Kubo relation for viscosity, as described in 
Sec. II, allows a calculation of viscosity using as its input 
a time series record of the motion of individual parti- 
cles jl3-15|. The relation is based on fluctuations, not a 
macroscopic velocity gradient. The Green-Kubo relation 
has been widely used in molecular dynamics (MD) sim- 
ulations, for example [l6j . However, there is a need for a 



test of the Grccn-Kubo relation for viscosity using an in- 
put of experimental data. In our search of the literature 
we have not found any such test, nor have we found any 
previous use of the Green-Kubo relation with an input of 
experimental data to determine viscosity. 

Viscosity is most commonly determined in experiments 
using a macroscopic velocity gradient. For example, in a 
rheometer 17], a stress is applied by a moving boundary, 
causing the liquid to flow with a velocity gradient, so that 
the viscosity can be determined using its hydrodynamical 
deflnition. Another experimental approach, which is used 
with experimental data from colloidal suspensions and 
other soft materials without macroscopic velocity gradi- 
ents, is the measurement of the mean-square displace- 
ment (MSD) of individual particles and the assumption 
that the Stokes-Einstein relation is valid, allowing a cal- 
culation of the viscosity from the measured MSD 3, 1^ , 
as discussed in Sec. II. All of these experimental methods 
are different from the use of the Green-Kubo relation. 

It has been questioned whether transport coefficients 
exist in 2D systems. Molecular dynamical simulations 
suggested that the diffusion, viscosity, and thermal con- 
ductivity coefficients would not exist in a 2D system of 
hard disks This result led to theoretical investiga- 

tions that indicated the time integral in the Green-Kubo 
relations diverges and the corresponding transport co- 
efficients in 2D systems are non-existing |2li. We will 
discuss this issue of convergence in Sec. IV E. 

In this paper, we test the Green-Kubo relation with 
an input of experimental data from a quasi-2D dusty 
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plasma. The dust particles in a dusty plasma, unlike 
hard disks, have a relatively long-range interparticle in- 
teraction. 

Our results and conclusions in this paper are intended 
to have a usefulness that extends beyond the area of 
dusty plasmas. Accordingly, we have attempted to make 
our presentation accessible to scientists who are not spe- 
cialists in that area. The experimental data used in this 
paper are from an experiment by Feng et al. We 
will introduce the concepts of dusty plasmas, along with 
some key aspects of the experiment, below as well as in 
Sec. III. Further details of the experiment can be found 
in @. 

A dusty plasma is a four- component mixture consist- 
ing of micron-size particles of solid matter, a gas of neu- 
tral atoms, free electrons, and free positive ions. We 
will refer in this paper to the particles of solid matter as 
"dust particles." In experiments, the gas is often argon, 
and so are the ions. The electrons and positive ions are 
present because of electric fields, provided by a power 
supply, which lead to a weak ionization of the neutral 
argon gas. Because of the much larger mobility of the 
electrons, as compared to the ions, many more electrons 
are collected on the dust particles, so that a dust particle 
develops a large negative electric charge, comparable to 
thousands of elementary charges [22[ . The electrons and 
ions that surround a negatively-charged dust particle are 
disturbed, resulting in a screening layer with a surplus of 
positive ions that surrounds the dust particle. Because 
of the complicated interactions amongst its four compo- 
nents, a dusty plasma is sometimes termed a "complex 
plasma" [i^l- While the gas, electrons and ions fill a 3D 
vacuum chamber, the dust particles do not. They are lev- 
itated against the downward force of gravity by a strong 
vertical electric field. In the experiment [8|, enough dust 
particles were introduced to fill a single horizontal layer, 
but not enough to begin filling a second layer. The dust 
particles were not in contact with any solid boundary, 
but they suffer a friction due to the surrounding neutral 
argon gas. In this paper, we will always consider the gas 
atoms as a whole to be a continuum, but we will consider 
the dust particles as individual entities. 

The experiment can be described as quasi-2D. The 
dust particles, although they are three dimensional, are 
arranged in two dimensions. It was verified, using video 
observations, that the dust particle motion in the verti- 
cal direction was extremely limited, as compared to that 
in the horizontal direction, so that dust particles moved 
past one another only as a result of their horizontal mo- 
tion. 

The neutral argon gas is rarefied in the experiment Q . 
At a density five orders of magnitude less than in a stan- 
dard atmosphere, collisions amongst gas atoms have a 
long mean free path, on the order of 1 cm. The effect of 
those collisions is diminished even more because the dust 
particles occupy only a thin 2D layer, so that a gas atom 



that collides with a dust particle is likely to be knocked 
out of the layer of the dust particles [2J|. Thus, no signif- 
icant transfer of momentum between two dust particles 
can occur due to the first one colliding with a gas atom 
which then collides with the second one. The only in- 
teraction between dust particles and gas that we must 
consider is the frictional drag force Ff on the dust parti- 
cle, which is proportional to the relative velocity between 
the dust particle and the g whole. In M the gas 

flow was negligibly slow. 

Since it is only the motion of the dust particles that 
will be of interest here, we will simplify our description 
of the four-component mixture. The dust particles are 
assumed to interact amongst themselves with a screened 
Coulomb repulsion, as discussed in Sec. IV. The role of 
the electrons and ions is then only to modify the inter- 
particle potential and provide the screening. Thus, in 
our simplified description of the four-component mixture 
we consider only a binary mixture: first, moving charged 
dust particles whose interaction potential is a screened 
Coulomb repulsion, and second a stationary neutral gas 
that exerts a frictional drag on moving dust particles. 

This reduction of a four-component mixture to a bi- 
nary mixture, in which all the properties of the electrons 
and ions are contained in the screening, has been used 
previously in theoretical descriptions of dusty plasmas, 
for example in the analysis of wave propagation [25| . 
When using this binary-mixture description, one could 
consider a charged dust particle as a "dressed parti- 
cle" [2^ consisting of a micron-size solid core that is neg- 
atively charged and a larger surrounding screening region 
that is positively charged. The center of this dressed par- 
ticle corresponds to what is observed experimentally by 
video microscopy. 

In the experiment [8j , the repulsion between dust par- 
ticles was so strong that the dust particles self-organized 
in a solid-like arrangement called a Wigner crystal [27 1. 
In order to study a liquid and its viscosity, this solid 
was melted by increasing the kinetic energy of the dust 
particles, which was done by using the laser-heating 
method 

When we refer in this paper to viscosity, it is only 
the motion of the dust particles that we directly take 
into account. In our simplified description of a dusty 
plasma, treating it as a binary mixture of dust particles 
and gas, we do not consider the momentum carried by 
electrons and ions. Moreover, a transfer of momentum 
between two dust particles does not occur due to the 
first dust particle colliding with a gas atom which then 
collides with the second, as discussed above. Thus, in 
our simplified binary-mixture description, the viscosity 
describes motion of only dust particles [29| . 

Previously, viscosity was studied in other dusty plasma 
experiments by applying a macroscopic shear stress us- 



ing laser manipulation |30l - l33| to generate a macroscopic 



velocity gradient, and using a hydrodynamical approach 
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to calculate the viscosity based on the measured velocity 
profile of the dust particles. In the test reported in this 
paper, we will compare the hydrodynamical result of [sij 
to the viscosity determined using a theoretical Green- 
Kubo relation with an input of data from the experiment 
of , which was performed without a macroscopic veloc- 
ity gradient. 

In Section II, the Green-Kubo relation for viscosity is 
reviewed. In Sec. Ill, we provide further details of the 
experiment Q. In Sec. IV, we introduce how we use the 
Green-Kubo relation with an input of experimental data. 
In Sec. V, we report our MD simulations of the exper- 
iment [ij. We determine the viscosity in Sec. VI, and 
in Sec. VII this result using the Green-Kubo relation is 
compared to the results of a previous experiment using 
a hydrodynamical method [3l|. In Sec. VII we also pro- 
vide a comparison to the results of previous computer 
simulations 2J, |3J, |35 1 . 



II. GREEN-KUBO RELATION 

To obtain transport coefficients such as diffusion, vis- 
cosity, and thermal conductivity in a liquid, Green-Kubo 
relations l^- 15 , [13, are often used. Their required 
inputs include time series for the positions and veloci- 
ties of particles. The simplest of these Green-Kubo re- 
lations is the one for diffusion. It can be derived easily 
using the physical assumption that the MSD for fluctu- 
ating particle displacements is proportional to the diffu- 
sion coefficient and the time |15|. The derivation of the 
Green-Kubo relation for viscosity is less trivial, a.nd it is 
based on the fluctuating stress, not an MSD [37|. Here 
we review the standard Green-Kubo relation for calcu- 
lating viscosity, in three steps, as it is used for all kinds 
of liquids, not just dusty plasmas in a liquid phase. 

In the first step, an off-diagonal element of the stress 
tensor Pxy{t) is defined by 
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where i and j denote different particles, N is the total 
number of particles of mass m, = {xi, yi) is the position 



of particle i. 



Xj, Uij — yi yj, 



and ^{rij) is the interparticle potential energy. Although 
not indicated in Eq. ([T]), the positions and velocities of 
particles vary with time, and this accounts for the time 
dependence of Pxy{t). In the second step, an autocorre- 
lation function of Pxy{t) is calculated as 



(2) 



We will refer to C,, (t) as the stress autocorrelation func- 
tion (SACF). The brackets (• • •) denote an average over 
an equilibrium ensemble, which in practice is often re- 
placed by an average over different initial conditions. In 



the third step, the SACF is integrated over time to yield 
the viscosity ry; for a 2D system the result is 
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Crj{t)dt, 



(3) 



where A is the area of the 2D system and T is its temper- 
ature. Equation ([3]) combined with Eq. ([2]) represent the 
Green-Kubo relation for viscosity in 2D. Similar Green- 
Kubo relations can be written for diffusion and thermal 
conductivity [15,, i36] . 

Viscosity rj and mass density p have different dimen- 
sions in 2D and 3D. The units of 77 are kgs~^ in 2D, and 
kgm~^s~^ in 3D. Correspondingly, in the denominator of 
Eq. ([3]) we have replaced the usual volume for a 3D sys- 
tem with an area A for the 2D system. In 2D, the areal 
mass density is p = nm with units of kgm^^, where n is 
the areal number density. We will report results for the 
kinematic viscosity rj/p, which has the same units in 2D 
and 3D. 

While Green-Kubo relations have been commonly used 
in computer simulations, their use with experimental 
data is uncommon. We are aware of only one previous 
calculation of any transport coefficient using the input of 
experimental data in a Green-Kubo relation. Using data 
from a dusty plasma experiment, Vaulina et al. [38l . Isoj 
obtained the diffusion coefficient using its Green-Kubo 
relation, which is a time integration of the velocity au- 
tocorrelation function. Calculating this autocorrelation 
function required an input of dust particle velocities, 
which were determined from experimental measurements 
of dust particle positions. In principle, the approach of 
Vaulina et al. of using a Green-Kubo relation to obtain 
the diffusion coefficient could be extended to other trans- 
port coefficients: the viscosity, thermal conductivity, and 
bulk viscosity. 

In this paper, we use the the Green-Kubo relation for 
viscosity with an input of experimental data. A compari- 
son of our result for viscosity with values determined in a 
previous experiment using a hydrodynamic method with 
a velocity gradient 31[ will serve as an experimental test 
of the Green-Kubo relation for viscosity. 

Besides the Green-Kubo relation, microrheology is an- 
other method to obtain the viscosity of fluids without 
a macroscopic velocity gradient [13 E^- In this ap- 
proach, the MSD of individual microparticles is mea- 
sured, and the Stokes-Einstein relation is assumed to 
be valid [H [3]. The Stokes-Einstein relation ^ is 
a combination of the Stokes law, which is a hydrody- 
namic model for viscous flow at a low Reynolds number, 
and the Einstein relation, which relates a diffusion coef- 
ficient for Brownian motion and a frictional force. This 
MSD-based method has been used in physical systems 
like colloidal suspensions, where a microparticles's mo- 
tion is overdamped due to the surrounding liquid sol- 
vent [Hill. 
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The Green-Kubo relation for viscosity is different from 
MSD-based approaches of determining viscosity. The 
derivation of the Green-Kubo relation centers on the fluc- 
tuations of the stress P^y, and it does not rely on the 
validity of a diffusion coefffcient. The Stokes-Einstein re- 
lation, while having great utility for many physical sys- 
tems, is known to fail for others such as supercooled liq- 
uids [41|]. Moreover, it is possible that a physical system 
can have a valid viscosity coefficient but lack a valid dif- 
fusion coefffcient, for example due to superdiffusion, as 
has been suggested for 2D systems, such as Yukawa liq- 
uids [illli]. 



III. EXPERIMENTAL INPUT 

Before reviewing the experiment Q we will discuss a 
few properties of dusty plasmas and their relevant length 
and time scales. 

When dust particles have a charge of several thousand 
elementary charges, their interparticle potential energy 
$ can be larger than their kinetic energy. In this case, 
the collection of dust particles is said to be a "strongly- 
coupled plasma" [42|, |43[ . A measure of strong coupling 
is the Coulomb coupling parameter F = (Q^ /AneoakBT), 
where Q is the charge of the dust particle, eo is the per- 
mittivity of free space, a is a typical interparticle distance 
as deffned below, and T is the kinetic temperature of the 
dust particles. The Coulomb coupling parameter is essen- 
tially a ratio of interparticle potential energy and kinetic 
energy. A plasma is strongly coupled when F > 1 , and it 
can behave like a liquid or a Wigner crystal. In a dusty 
plasma, the dust particles are usually strongly coupled 
due to their large charge Q. The electrons and the ions 
have a much smaller charge, and in most plasmas they 
are not strongly coupled, unless great efforts are made to 
cool them to low temperatures [J|. 

The length scales of a dusty plasma include the screen- 
ing length Xjj and the typical distance between dust par- 
ticles. In the dusty plasma literature, the typical distance 
between dust particles is commonly reported either as the 
lattice constant b for a Wigner crystal, or as a = (nn)^^/^ 
for a liquid, where n is the areal number density. In the 
literature for strongly-coupled plasmas, a is called the 
2D Wigner-Seitz radius [4^. All three of these length 
scales, Xd, b and a, are typically on the order of 1 mm. 
The diameter of a dust particle is typically a few mi- 
crons, which is much smaller than any of these length 
scales, and also much smaller than the mean free path 
for collisions amongst the rarefied argon gas atoms. 

The time scales of a dusty plasma include measures 
of collective motion amongst the dust particles, and 
of the frictional drag experienced by a dust particle 
due to the g whole. The former is quantified 

by the 2D plasma frequency for dust particles, ujpd = 
{Q^ /2neQma'^y/^ [42], where the subscripts p and d refer 



to plasma and dust, respectively. The frictional drag 
experienced by a dust particle due to the neutral ar- 
gon gas is quantified by a friction coefficient, which is 
commonly defined in the literature for dusty plasmas as 
h'f = Ff/mv, where Fj is the gas friction force experi- 
enced by one dust particle and mv is the momentum of 
the same particle. This coefficient has the dimension of 
inverse time, and a typical value for a micron-size dust 
particle in a rarefied gas is z^/ « 1 s^^. If the gas is 
sufficiently rarefied, ujpd > Vf, so that the dust particle 
motion is said to be underdamped. 

We now summarize some aspects of the experiment Q 
that are relevant for our analysis. The dust particles 
were polymer microspheres of 8.1 /im diameter. They 
were levitated by a vertical dc electric field to form a 
single horizontal layer, as sketched in Fig. 1. Radial dc 
electric fields provided horizontal confinement so that the 
dust particles filled a circular region of diameter 52 mm. 

For each experimental run, one movie of dust parti- 
cle motion was recorded. A Phantom v5.2 high-speed 
camera viewed the dust particles from above. It was op- 
erated at 250 frames/s, so that data were recorded with 
a time interval At = 4 ms. The duration of a movie, 
20 s, was limited by the camera's memory. The cam- 
era's field-of-view (FOV) included N « 2100 dust parti- 
cles. Two of the required inputs for Eq. ([T]), the positions 
and velocities of the dust particles, were obtained using 
the moment method 44 1 and particle tracking [iH . The 
interparticle distance, averaged over the camera's FOV, 
was characterized by 6 = 0.67 mm in the Wigner crystal, 
and a — 0.35 mm in the liquid. 

A total of seven runs were performed. Three runs with- 
out laser heating were made to determine the interpar- 
ticle potential energy in the Wigner crystal. Four runs 
were done with laser heating to make the liquid, and they 
provide the data that we will use here as the input for 
the Green-Kubo relation, to determine viscosity in the 
liquid phase. 

The laser heating method [3] uses the radiation- 
pressure due to laser beams, which are directed toward 
the dust particles by scanning mirrors 47[, which we 
configured as in [23, li^. Laser heating increases the 
kinetic temperature of dust particles without changing 
their charge. The velocity distribution function has been 
observed to be nearly Maxwellian [28|. In the experi- 
ment Q, the kinetic t emp erature T, calculated from the 
mean-square velocity [49|, was 2.5 x 10* K in the hquid 
(with laser heating), and 10^ K in the Wigner crystal 
(without laser heating) . The temperature was nearly uni- 
form spatially, which is desirable because viscosity varies 
with temperature. 
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IV. USING THE GREEN-KUBO RELATION 
WITH EXPERIMENTAL INPUT DATA 

Calculations of viscosity with the Green-Kubo relation 
Eqs. (P)-® must be adapted to use the experimental 
data as input, due to four difficulties. First, the experi- 
ment Q provided direct measurements of Xi, yi, Vix, and 
Viy, but not of the interparticle potential energy ^(ry ). 
Second, the camera has a finite FOV, so that we do not 
have data for all dust particles in the 2D layer. Third, the 
motion of dust particles may include a local macroscopic 
flow, i.e., a non-zero time- average velocity. Fourth, the 
data for the positions and velocities of the dust parti- 
cles are recorded as a time series of a finite duration, so 
that the integral in Eq. must have a finite limit. We 
will take all this into account by making a number of 
approximations in the Green-Kubo calculations. We will 
next describe these approximations as well as discuss the 
validity of the results that are obtained. 



A. Interparticle potential 

To solve the first difficulty, a lack of experimental mea- 
surements of the interparticle potential energy ^(rij), 
we will use a model for these energies when calculat- 
ing the off-diagonal element of the stress tensor P^y For 
a single-layer dusty plasma like ours, models that have 
been tested successfully include an isotropic repulsion ac- 
cording to the 3D Debye-Hiickel potential [50] 



Q exp(-r^j/Ag) 
47reo Ti, 



(4) 



as well as more complicated isotropic [5l| and non- 
isotropic interactions [5^. Here we will use the Debye- 
Hiickel potential, Eq. (jlj, where i and j are dust particles 
of charge Q separated by a distance r^, and A_d is the 
screening length due to electrons and ions. The corre- 
sponding potential energy is = Q(f>{rij). In the lit- 
erature for dusty plasmas, it is common to name Eq. (|4]) 
after Yukawa rather than Debye and Hiickel. This po- 
tential has been used in theoretical and simulation stud- 
ies of visc osity in strongly coupled plasmas, for exam- 
ple [iilsi-IHei. 

Two parameters in Eq. (HI), Q and Xd, were de- 
termined in the experiment [8| using a phonon-spectra 



method [57|, |58[ in the Wigner crystal. For the three ex- 
perimental runs without laser heating, position and ve- 
locity measurements were used to compute the phonon 
spectra, which were compared to theoretical wave dis- 
persion relations [2^ that assume Eq. Q , yielding Xd = 
0.70 ± 0.14 mm and Q/e = -6000 ± 600, where e is the 
elementary charge. Using these two values along with 
measured particle positions and velocities, we can calcu- 
late Pxy and then the SACF j59i] . 



Having determined Q and Ad, we can calculate the 
values of other parameters. We find that the Coulomb 
coupling parameter is F = 68 for the liquid of dust par- 



ticles. The dimensionless particle spacing [42| is k = 
a/Xo — 0.5 ± 0.1. We also calculated ujpd = 30 ± 3 s^^. 
We note that ujpd is significantly larger than the gas fric- 
tion coefficient ~ 2.4 s^^, indicating that dust particle 
motion is undcrdamped. 



B. Finite field of view 

To solve our second difficulty, the finite field of view 
(FOV) of the camera that limits us to recording data for 
only a portion of the dust particles, we will cut off the 
interparticle potential in Eq. (|4]) at large distances. In 
addition, we will divide the FOV into inner and outer 
regions. Fig. 2. The cutoff is done at a large interparticle 
distance of 56, where the exponential in Eq. ([4]) is < 10^^. 

The FOV is divided into inner and outer regions be- 
cause the potential energy of a dust particle cannot be 
obtained meaningfully if it is located near the edge of the 
FOV, due to interactions with dust particles of unknown 
positions outside the FOV. Therefore, dust particles in 
the outer region are used only to calculate the potential 
energies $ of dust particles in the inner region. In other 
words, when calculating P^y, we limit the dust particles 
i to those that are located within the inner region, and 
we account for their interaction with other dust particles 
j located in both the inner and outer regions, as shown in 
Fig. 2. The outer region has a width of 56 to allow for a 
cutoff radius of 56, and the inner region is 22.06 x 23.76 
(i.e., 14.8 X 15.9 mm^). Because of this different treat- 
ment of dust particles in the inner and outer regions, we 
calculate P^y as 



Pxy — 



^ix ^xy 



1 ^ 



(5) 



where the stress Pxy is (implicitly) a function of time. 
All time series data are recorded at a time interval At = 
4 ms. For the experiment M « 600 dust particles are 
in the inner region, while N « 2100 are in both regions 
combined. Both N and A/ fluctuate slightly with time, 
as dust particles move across the edges of the regions, but 
not enough to affect the result for viscosity significantly. 



C. Non-zero time averages 

Our third difficulty to solve is that the average velocity 
is not zero, due to finite macroscopic fiow velocities of the 
dust particles, despite efforts that were made in the ex- 
periment Q to avoid them. When computing the stress, 
the dust particle velocities are assumed to fluctuate about 
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an average value of zero. In fact, a non-zero average ve- 
locity would contribute an unwanted constant value to 
the Pxy time series, when computed using Eq. ([!]) or ([5]), 
which would cause the SACF, C^, to decay to a non-zero 
value and introduce an unphysical contribution to the 
calculated viscosity. To solve this difficulty, without any 
approximation, we subtract from Pxyit) its time-average 
value, yielding then the fluctuating portion Pxy{t)- We 
then replace Eq. ^ for the SACF by 



^vi^) ~ {Pxy{t + to)Pxy{to))o- 



(6) 



Here, the brackets (• • •)o in Eq. ([6]) denote an average 
over various initial times to, if the data are from a single 
run. 

Our results for the SACF, for the four runs with laser 
heating, are shown in Fig. 3. All four runs show the 
same general trends, which resemble those seen in MD 
simulations as in 34| . 



The collection of dust particles is best de- 

48j . where the 



D. Integration limit 



equilibrium 

scribed as a driven-dissipative system 
driving is mainly provided by the laser beams, and the 
dissipation is provided by gas-dust collisions as well as 
the dust-particle viscosity. Despite these nonequilibrium 
conditions, however, the velocity distribution function 
for the dust particles has been observed to be nearly 
Maxwellian [28|. Thus, the statistics are close to those 
of an equilibrium system, which motivates us to use the 
Green-Kubo relation. 

Third, we must ask whether long-time tails in the cor- 
relation function prevent the convergence of its integral, 
as was predicted theoretically 2l| for 2D systems with 
hard-disk interparticle interactions. More recently, for 
a 2D system with a Debye-Hiickel potential, simulations 
by Donko et al. 36| indicated that long-time tails occur 
in some but not all cases. In particular, they reported 
that the SACF decays fast enough that its time integral 
77 converges for a liquid at temperatures near the melting 
point, but not at absolute temperatures far above the 
melting point. 



To solve the fourth difficulty, the finite time duration 
of data, the viscosity 77 is computed with a finite upper 
limit in the time integral of the SACF. In principle the 
upper limit should be infinite, as in Eq. ([3|), but we use 
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AkeT 



(7) 



and we follow the practice used in MD simulations of 
choosing the integration limit i/ as the time when C,,(t) 
crosses zero [H^, as shown in Fig. 3. Because our time 
series C,,(t) is noisy, we count a zero crossing only if it 
results in C^(t) remaining negative for at least 2/S.t. In 
Eq. ([7]) , v4 is the area of the inner region. 



E. Validity 

For using the Green-Kubo relation to calculate viscos- 
ity with experimental input data, we should ask whether 
the approach is valid, when it is used for a dusty plasma. 
We will mention three questions. 

First, we note that the Green-Kubo relations are 
strictly speaking for the thermodynamic limit, where the 
number of dust particles and the system size tend to 
infinity while keeping the number density constant. In 
fact, the experiment has only thousands of dust parti- 
cles. However, we believe that our experimental system 
size is large enough to use a Green-Kubo relation, as in- 
dicated by our system-size tests in Sec. V. 

Second, we must consider the distinction between equi- 
librium and nonequilibrium systems. While our laser- 
heated dusty plasma is in a steady state, it is not in 



V. SIMULATIONS 

To assess three sources of error mentioned in Sec. IV, 
viz. the potential cutoff, the FOV division, and the fi- 
nite system size, we used MD simulations based on the 
Langevin equation 61|. In these simulations we numeri- 
cally integrate the Langevin equation to obtain the mo- 
tion of each dust particle. This equation gives the force 
acting on a dust particle as three terms: a sum of an elec- 
tric force due to all other dust particles using Eq. (jH), 
a mean friction Ff due to the gas as a whole, as well 
as Gaussian random forces around this mean, to model 
the collisions of the dust particle with the neutral gas 

63, m. 



atoms 



Dust particle positions, velocities, and interaction en- 
ergies were recorded at each time step of 0.019 . We 
used N — 4096 particles in a 2D rectangular box with 
periodic boundary conditions. In a 2D Langevin MD 
simulation with a Debye-Hiickel potential like ours, the 
equations of motion have three dimensionless parameters: 
the friction coefficient Vf /ujpd, the Coulomb coupling pa- 
rameter r, and the dimensionless particle spacing k. To 
mimic the experiment [8|, we used Uf /ujpd = 0.08, V = 68, 
and K = 0.5. 

To test the effect of the cutoff, we carried out simula- 
tions at two cutoff distances to estimate the systematic 
error introduced in the last term of Eq. ([SJ. We found 
that for a cutoff of 56, which we use in this paper, the vis- 
cosity result rj was reduced by less than 5%, as compared 
to the result for a much larger cutoff of 136. 

To test for errors arising from the division of the FOV, 
we compared results with and without the division, using 
Eq. ([S]) and Eq. ([ij, respectively. We found the viscosity 
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differed only negligibly. 

Finally, to test for the effects of finite system size, we 
compared our Langevin MD simulation results for two 
sizes. Results for the larger system size of 4096 particles 
are reported in 2J]. For the smaller size, to mimic the 
experiment [§| , we used 48 simulation runs with M = 600 
particles, as in the inner region for the experiment 0, for 
a duration [6^ of 607 . Comparing these two results 
for 77, we find no statistically-significant system-size ef- 
fects. This test, shown in Table I, gives us confidence 
that the number of dust particles in [8| was not so small 
as to preclude using the Green-Kubo relation. 



was much stronger than gas friction, Wpd ^ r^/, in both 
experiments. A difference in the experimental conditions 
was that the two laser beams were manipulated differ- 
ently so that in [3l| they produced a macroscopic ve- 
locity gradient, while in [8| there was heating without a 
macroscopic velocity gradient. 

Considering the complexity of the dusty plasma and 
our simplifications in describing it. we cannot expect our 
result to agree with the results of [31| better than within 
about a factor of two. In fact, we find agreement within 
this factor when comparing our result rj/ p = 0.16 ± 0.02 
using the Green-Kubo relation and the range of values 
reported by Nosenko and Goree [sjj. Fig. 4. 



VI. RESULT FOR VISCOSITY USING THE 
GREEN-KUBO RELATION WITH 
EXPERIMENTAL DATA 

We now present our result for the viscosity 77 using 
the Green-Kubo relation, Eqs. ©-(III). We report the 
kinematic viscosity, 7]/ p, to allow convenient comparison 
to other experiments and simulations. 

We find -q/p = 0.16 ± 0.02, in units of a'^ujpd- The 
value of the normalization factor for the experiment |;8] 
is a^ujpd — 3.7 X 10~® m^/s, while the areal mass den- 
sity is p = 1.1 X 10^^ kg/m^. The value of 0.16 is the 
mean of the four experimental runs in the presence of 
laser heating, as plotted in Fig. 3. The error estimate of 
± 0.02, calculated as the standard deviation of the mean, 
indicates the run-to-run random variation. 



VII. DISCUSSION 
A. Test of the Green-Kubo relation 

Our most significant result is a test of the Green- 
Kubo relation for viscosity, using an input of experimen- 
tal data. We perform this test by comparing our result 
J]/ p = 0.16 ± 0.02 to the previously-reported result, from 
a dusty plasma experiment fsif that used a hydrodynam- 
ical approach. 



In the experiment of Nosenko and Goree [3l| , the dust 
particles fiowed in their quasi-2D layer with a macro- 
scopic velocity gradient, allowing a determination of vis- 
cosity using a Navier Stokes equation of motion for the 
local flow velocity of the dust particles. Differently 
from [31I, in our analysis here we consider the dust par- 
ticles only as individual particles, while in [3l| data for 
individual particles were averaged to allow considering 
the dust particles as a continuum, as is necessary for 
a hydrodynamical approach. Both experiments, [s!] and 
[sif . were performed in the same chamber, and both had 
a quasi-2D layer of dust particles. The reported values 
of K were nearly the same, k = 0.53 for 31| and 0.5 
for The electrical interactions among dust particles 



B. Comparison to simulations 

We also compare our result from the Green-Kubo re- 
lation to the available data for viscosity from the simu- 
lation literature. Fig. 5. All these data are from 2D MD 
simulations with a Debye-Hiickel potential, and most of 
them use the Green-Kubo relation, except the nonequi- 
librium simulation of [35| , which used a so-called nonequi- 
librium simulation method to produce a macroscopic ve- 
locity gradient. We find that the simulations predicted 
values that are larger, by about a factor of two, than our 
result for the Green-Kubo relation using experimental 
input. 

The discrepancy between our result here and the sim- 
ulation results in Fig. 5 could arise from the differences 
between the simulations and the experiment These 
differences include the use of periodic boundary condi- 
tions in the simulations to mimic an infinite system, while 
the dust particles in the experiment fill a finite region due 
to dc radial electric fields. While all particles in the sim- 
ulations are identical, those in the experiment vary in 
diameter by a few percent [65| . with a similar variation 
in charge. Nonuniformities occur in the simulations only 
as transients due to fluctuations, while the experiment 
has static nonuniformities due to the dc radial electric 
fields. These dc fields induce static stresses that can lead 
to more structural disorder which can result in an easier 
deformation of the arrangement of dust particles and a 
reduction in the stress required to generate shear fiow, 
which is equivalent to a reduction of the viscosity. 

Another possible explanation for the discrepancy with 
the simulation results is that the potential in the exper- 
iment [§| may not be exactly a Debye-Hiickel potential, 
as assumed in Sec. IV A. One alternative, instead of as- 
suming a particular form for the potential, could be an 
empirically determined potential of mean force, calcu- 
lated from an experimentally-measured pair correlation 
function, as has been proposed theoretically 3, 66|. An 
advantage of this approach is that all physical processes 
that affect the potential would be included in the empir- 
ical result [66j]. 
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VIII. CONCLUSIONS 

The Green-Kubo relation for viscosity has been tested 
using an input of experimental data. The value for the 
viscosity determined by the Green-Kubo relation with 
the input of data from an experiment [8J was compared 
to the value from aprevious experiment using a hydrody- 
namical method [31]. In both experiments the physical 
system was a quasi-2D dusty plasma, and the conditions 
were similar, aside from the absence of a macroscopic ve- 
locity gradient in the experiment 3 for the Green-Kubo 
result. We found that the results agree as well as can 
be expected. This agreement serves as our test of the 
Green-Kubo relation for viscosity of a dusty plasma. 

Additionally, we compared our result for the viscosity 
determined by the Green-Kubo relation with the predic- 
tions of MD simulations 24, yj, bsl. The results were 



as consistent as expected, given the differences between 
the simulations and the experiment Q that provided our 
input data. 

Further tests are needed for other Green-Kubo rela- 
tions. Because the Green-Kubo relations for the various 
transport coefficients are all different, a test of the Green- 
Kubo relation for viscosity, as we have presented here, 
does not also serve as a test of another Green-Kubo rela- 
tion. In addition, tests for other physical systems, such 
as the 2D systems mentioned in Sec. I, would be useful. 

This work was supported by NSF and NASA. 



[1] 
[2] 
[3] 
[4] 
[5] 

[6] 

[7] 
[8] 

[9] 

[10] 

[11] 
[12] 
[13] 

[14] 
[15] 



Electronic address: yan-feng@uiowa.edu| 

C.A. Murray, W.O. Sprenger, and R.A. Wenk, Phys. Rev. 

B 42, 688 (1990). 

P.M. Reis, R.A. Ingale, and M.D. Shattuck, Phys. Rev. 
Lett. 96, 258001 (2006). 

C.C. Grimes and G. Adams, Phys. Rev. Lett. 42, 
795 (1979). 

T.B. Mitchell, J.J. Bollinger, X.-P. Huang, W.M. Itano, 
and D.H.E. Dubin, Phys. Plasmas 6, 1751 (1999). 
P.K. Shukla and A. A. Mamun, Introduction to Dusty 
Plasma Physics (Institute of Physics Publishing, Bristol 
and Philadelphia, 2002). 

G.E. MorfiU and A.V. Ivlev, Rev. Mod. Phys. 81, 
1353 (2009). 

A. Piel, Plasma Physics (Springer, Heidelberg, 2010). 
Y. Feng, J. Goree, and B. Liu, Phys. Rev. Lett. 105, 
025002 (2010). 

G.K. Batchelor, An Introduction to Fluid Dynamics 

(Cambridge University Press, Cambridge, 1967). 

M.S. Murillo, Phys. Rev. Lett 85, 2514 (2000). 

G. Salin, Phys. Plasmas 14, 082316 (2007). 

P.K. Kaw and A. Sen, Phys. Plasmas 5, 3552 (1998). 

M.S. Green, J. Chem. Phys. 20, 1281 (1952); 22, 398 

(1954). 

R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 

J. P. Hansen and I.R. McDonald, The Theory of Simple 

Liquids, 2nd ed. (Elsevier Academic Press, Amsterdam, 



[16] 

[17 
[18 

[19 



[20 
[21 

[22' 
[2i 

[24 

[25; 

[26 

[27 

[28; 

[29 



[30] 
[31] 
[32] 
[33] 

[34] 
[35] 

[36] 

[37] 



1986). 

Z. Donko, P. Hartmann, and J. Goree, Mod. Phys. Lett. 
B 21, 1357 (2007). 

R. Pal, J. Colloid Interface Sci. 225, 359 (2000). 

T.G. Mason and D.A. Weitz, Phys. Rev. Lett. 74, 1250 

(1995). 

J.C. Crocker, M.T. Valentine, E.R. Weeks, T. Gisler, 
P.D. Kaplan, A.G. Yodh, D.A. Weitz, Phys. Rev. Lett. 

85, 888 (2000). 

B. J. Alder and T.E. Wainwright, Phys. Rev. Lett. 18, 
988 (1967), and Phys. Rev. A 1, 18 (1970). 

M.H. Ernst, E.H. Hauge, and J.M.J, van Leeuwen, Phys. 
Rev. Lett. 25, 1254 (1970). J.R. Dorfman and E.G.D. Co- 
hen, Phys. Rev. Lett. 25, 1257 (1970). 

C. K. Goertz, Rev. Geophys. 27, 271 (1989). 

D. Samsonov, J. Goree, H.M. Thomas, and G.E. MorfiU, 
Phys. Rev. E 61, 5557 (2000). 

Y. Feng, J. Goree, and B. Liu, Phys. Plasmas 18, 057301 
(2011). 

X. Wang, A. Bhattacharjee, and S. Hu, Phys. Rev. Lett. 

86, 2569 (2001). 

G. Joyce, M. Lampe, and G. Ganguh, IEEE Trans. 
Plasma Sci. 29, 238 (2001). 

L. Bonsall and A.A. Maradudin, Phys. Rev. B 15, 1959 
(1977). 

V. Nosenko, J. Goree, and A. Piel, Phys. Plasmas 13, 
032106 (2006). 

The dilTerent roles of the dust particles and the neutral 
gas in determining the transfer of momentum can be illus- 
trated by comparing a dusty plasma to a dilute colloidal 
suspension. The latter consists of two components: solid 
particles and a liquid solvent such as water, which have 
comparable mass densities. If the solid particles are suf- 
ficiently small, they are dragged along in the fiow of the 
solvent, so that the local velocities of the solid particles 
and solvent are nearly the same, G.K. Batchelor, op. cit. 
Then, the two components of a dilute colloidal suspen- 
sion can be considered to be moving as one, which has an 
effective viscosity that is slightly larger than that of the 
pure solvent according to Einstein's result: A. Einstein, 
Ann. Phys. 34, 591 (1911). A dusty plasma is different 
because instead of a liquid solvent there is a rarefied gas, 
which is many orders of magnitude less dense. Then the 
dust particles can fiow with a velocity that is different 
from that of the rarefied gas as a whole. This difference 
in velocity between a dust particle and the gas as a whole 
leads to a friction between them, which is different from 
the viscous friction due to collisions amongst the dust 
particles via their electrical repulsion. 
W.-T. Juan, M.-H. Chen, and L. I, Phys. Rev. Lett. 64, 
016402 (2001). 
V. Nosenko and J. 
(2004). 

A. Gavrikov, I. Shakhova, A. Ivanov, O. Petrov, 
N. Vorona, and V. Fortov, Phys. Lett. A 336, 378 (2005). 
N.A. Vorona, A.V. Gavrikov, A.S. Ivanov, O.F. Petrov, 
V.E. Fortov, and LA. Shakhova, J. Exp. Theor. Phys. 
105, 824 (2007). 

B. Liu and J. Goree, Phys. Rev. Lett. 94, 185002 (2005). 
Z. Donko, J. Goree, P. Hartmann, and K. Kutasi, Phys. 
Rev. Lett. 96, 145003 (2006). 

Z. Donko, J. Goree, P. Hartmann, and B. Liu, Phys. Rev. 
E 79, 026401 (2009). 

D.A. McQuarrie, Statistical Mechanics. (Harper & Row, 



Goree, Phys. Rev. Lett. 93, 155004 



9 



New York, 1976). 

[38] O.S. Vaulina, X.G. Adamovich, O.F. Pctrov, and 
V.E. Fortov, Phys. Rev. E 77, 066404 (2008). 

[39] O.S. Vaulina, X.G. Adamovich, O.F. Petrov, and 
V.E. Fortov, J. Exp. Theor. Phys. 107, 313 (2008). 

[40] R.A.L. Jones, Soft Condensed Matter. (Oxford University 
Press, Auckland, 2003). 

[41] E. Rosslcr, Phys. Rev. Lett. 65, 1595 (1990). 

[42] G.J. Kalman, P. Hartmann, Z. Donko, and M. Rosen- 
berg, Phys. Rev. Lett. 92, 065001 (2004). 

[43] S. Ichimaru, Rev. Mod. Phys. 54, 1017 (1982). 

[44] Y. Feng, J. Goree, and B. Liu, Rev. Sci. Instrum. 78, 
053704 (2007). 

[45] Y. Feng, J. Goree, and B. Liu, Rev. Sci. Instrum. 82, 
053707 (2011). 

[46] M. Wolter and A. Melzer, Phys. Rev. E 71, 036414 
(2005). 

[47] A. Melzer, S. Nunomura, D. Samsonov. Z.W. Ma, and 

J. Goree, Phys. Rev. E 62, 4162 (2000). 
[48] B. Liu and J. Goree, Phys. Rev. Lett. 100, 055003 (2008). 
[49] Y. Feng, J. Goree, and B. Liu, Phys. Rev. Lett. 100, 

205007 (2008). 

[50] U. Konopka, G. E. Morfill, and L. Ratke, Phys. Rev. 
Lett. 84, 891 (2000). 

[51] O. S. Vaulina, E. A. Lisin, A. V. Gavrikov, O. F. Petrov, 
and V. E. Fortov, JETP, 110, 662 (2010). 

[52] Non-isotropic models have been developed to take into 
account the effects of ion wakes, which arise from a 
flow of ions that is driven by the same vertical electric 
field that levitates the dust particles. Examples include 
V.A. Schweigert, I.V. Schweigert, A. Melzer, A. Homann, 
and A. Piel, Phys. Rev. Lett. 80, 5345 (1998); M. Lampe, 
G. Joyce, G. Ganguli, and V. Gavrishchaka, Phys. Plas- 



mas 7, 3851 (2000); and R. Kompanects, U. Konopka, 
A.V. Ivlcv, U. Tsytovich, and G. Morfill, Phys. Plasmas 
14, 052108 (2007). 

[53] M.S. Murillo, Phys. Rev. E 62, 4115 (2000). 

[54] T. Saigo and S. Hamaguchi, Phys. Plasmas 9, 1210 

(2002) . 

[55] G. Salin and J.-M. Caillol, Phys. Plasmas 10, 1220 

(2003) . 

[56] G. Faussurier, Phys. Rev. E 69, 066402 (2004). 

[57] S. Nunomura, J. Goree, S. Hu, X. Wang, A. Bhattachar- 

jee, and K. Avinash, Phys. Rev. Lett. 89, 035001 (2002). 
[58] S. Nunomura, S. Zhdanov, D. Samsonov, and G. Morfill, 

Phys. Rev. Lett. 94, 045001 (2005). 
[59] We carried out tests that showed that the SACF is more 

sensitive to an error in Q than in Ad. 
[60] Z. Donko, J. Goree, and P. Hartmann, Phys. Rev. E 81, 

056404 (2010). 

[61] F. Reif, Fundamentals of Statistical and Thermal 
Physics. (McGraw Hill, New York, 1965), Chap. 15. 

[62] W. F. van Gunsteren and H. J. C. Berendsen, Molecular 
Physics 45, 637 (1982). 

[63] T. Ott and M. Bonitz, Phys. Rev. Lett. 103, 195001 
(2009). 

[64] For comparison, the memory time of collective motion in 

this liquid of dust particles was « 10 (^~^ , as measured 
by the zero-crossing of the correlation function Cr){t). 

[65] B. Liu, J. Goree, V. Nosenko, and L. Boufendi, Phys. 
Plasmas 10, 9 (2003). 

[66] U. de Angelis, Phys. Plasmas 8, 1751 (2001). 



10 



TABLE L Results for the kinematic viscosity 7;/p, which is normalized here by a^ujpd to make it dimensionless. Viscosity rj is 
reported as the mean for multiple runs; the standard deviation a and standard deviation of the mean cm are listed. 



data source 


calculation 




data size 


ri/p, units {a^copd) 




procedure 


M 


run duration 


runs mean a um 



experiment [8] Eq.(5)-i'(6)^(7) ^ 600 ?a 2100 ^ 607aj~j^ (20.2 sec) 4 0.16 0.04 0.02 

Langevin simulation (Sec. V) Eq.(5)-^(6H(7) 600 4096 607tJ-^^ 48 0.24 0.08 0.01 

equilibrium simulation [24] Eq.(l)^(6)-)-(7) 4096 4096 1.86 x lO^cj"^^ 4 0.26 0.02 0.01 

Langevin simulation [24] Eq.(l)-s-(6)^(7) 4096 4096 1.86 x lO^cj"/ 4 0.27 0.02 0.01 



FIG. 1; (color online) Configuration for the experiment [^. 
(a) This diagram of the vacuum chamber is shown in an ex- 
ploded view, to better show the lower electrode, which in the 
experiment was located inside the chamber. A low-pressure 
gas of neutral argon atoms filled the chamber. A radiofre- 
quency voltage was applied between two electrodes, sepa- 
rated by an insulator. One electrode was the powered lower 
electrode, and the other consisted of the grounded vacuum 
chamber and shield. The gas was partially ionized, yielding a 
plasma with three components: electrons, positive argon ions, 
and neutral argon atoms. The x and y axes correspond to the 
two orthogonal directions used in measuring the positions and 
velocities of dust particles. The side ports were used to admit 
laser beams, not shown here. For heating, a pair of 532 nm 
laser beams were directed into the chamber by scanning mir- 
rors as in [2^ , while for illumination a 577 nm laser sheet was 
used with a configuration as in [43|. (b) This sketch shows a 
side view of the lower electrode. Polymer microspheres were 
introduced by shaking them into the plasma from above, and 
they moved downward due to gravity g. They gained a neg- 
ative electric charge Q and were levitated upward due to a 
vertical dc electric field E so that they remained in a single 
horizontal layer above the lower electrode. There was also 
a smaller radial dc electric field, not shown, which provided 
horizontal confinement. Images of the dust particles were 
recorded by a video camera viewing through the top port. 
Shown here are 5 mm x 5 mm portions of two images: (c) a 
Wigner lattice without laser heating (lines have been drawn 
to indicate the lattice structure) and (d) a liquid maintained 
by laser heating. 



FIG. 2: (color online). Sketch of the division of the camera 
FOV into inner and outer regions for the experiment §]. In 
Eq. ((5]), the subscripts i and j refer to particles located in 
the inner region, and both the inner and outer regions, re- 
spectively. The circle indicates the cutoff distance 5& for the 
potential. The unused portion of the camera FOV on the 
right is not included in the analysis. 
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FIG. 3: (color online). Stress autocorrelation function SACF, 
i.e., Crfit), and its time integration. Results are shown for 
four runs of the experiment [3 . All quantities shown are nor- 
malized. The SACF is normalized as C = Cn/{AkBTpa^ujpj) 
and is drawn here at 10 x magnification. Time is given in 
units of oj~^, that is, r = tujpd. The integral shown by the 
smooth curve is used in the Green-Kubo relation, Eq. ([T]), to 
calculate the viscosity. Choosing the integration limit ti as 
the time when C(r) crosses zero as described in Sec. IV D, 
the integral yields the dimensionless viscosity, as indicated by 
the solid circle for each run. 



FIG. 4: (color online). Our value of the kinematic viscosity 
calculated from the Green-Kubo relation, using input from 
the experiment 8], compared to results from a hydrodynam- 
ical analysis of a previous experiment fsil]. Values are made 
dimensionless by normalizing by a^Upd- The x axis, which 
has a logarithmic scale, is the Coulomb coupling parameter F 
as defined in Sec. III. Our result, shown as a solid diamond, is 
the mean for four experimental runs in Fig. 3, and the vertical 
error bar indicates only the run-to-run variation, calculated 
as the standard deviation of the mean. The horizontal error 
bar (for the result from the Green-Kubo relation) reflects a 
10% uncertainty in Q. 



FIG. 5: (color online). Comparison of our result using the 
Green-Kubo relation for viscosity with the input of data from 
the experiment 0], shown as a diamond as in Fig. 4, to values 
from previously-reported 2D Debye-Hiickel simulations [2^ . 
[33. 135I]. The simulation results shown for [2^ are also listed 
in Table I. Both axes are logarithmic. 
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